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We introduce an unfamiliar Kondo ladder problem by coupling the well-known Takhtajan-Babujan 
S=l chain to a half-filled one-dimensional electron gas, and we solve it using quantum field theory 
techniques. Based on the so-called underscreened Kondo effect, we find an insulating spin liquid 
with both deconfined spinons and optical magnons. Important consequences are discussed. 

PACS numbers; 71.iO.Pm, 71.27.+a, 75.10.Jm 
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The problem of the co-existence between conduction 
electrons and an array of local magnetic moments re- 
mains an eminent unsolved problem in condensed matter 
physics, in link with the well-known Kondo effect. This 
poses a real challenge to the theorist. As a first step, the 
one-dimensional (ID) model of a Kondo chain with an 
M-fold degenerate band of conduction electrons interact- 
ing with a spin S Heisenberg chain has been solved in a 
semi-classical limit |l|] . At low energy, the system is insu- 
lating and the spin properties can be described by an 0(3) 
nonlinear sigma model (NLcrM) with a topological term 
e = tt{2S ~ M). For \M - 2S\ = (even), single electron 
excitations at low energies are massive spin polarons. For 
|M — 2S'| — (odd), the NLcrM in contrast becomes criti- 
cal and belongs to the universality class of the isotropic 
S=l/2 Heisenberg chain. In this Letter, we investigate 
the interesting case S=l and M=l in a so-called 'ladder' 
picture |g]. As we shall see below, the characterization 
of the spin array as an S=l quantum spin chain with 
prevalent gapless triplet excitations should allow to push 
forward our understanding of the underscreened Kondo 
effect H in the Kondo lattice model (KLM). For exper- 
imental applicability, this model appears to be partic- 
ularly appropriate to describe heavy-fermion materials 
with Uranium compounds such as UPt^ (as pointed out 
in ref. M, where we have firstly solved the problem of 
two localized S=l spins embedded in a host metal). 

Let us start with the integrable Takhtajan-Babujan 
S=l chain on a lattice, described by: Ti = 
■^l|E,(S,S,+i) - /3(S,S,+02 ^here /3 = 1 |. S, de- 
scribe local moments. The unusual term quartic in spin 
can be generated, for example, by phonons ||6|. Solv- 
ing Bethe Ansatz equations, the only elementary excita- 
tion is known to be a doublet of gapless spin- 1/2 spin 
waves, with total spin or 1 [g|. The approach dis- 
cussed in this Letter is based on field theory techniques 
like the non-Abelian Q and Abelian [|| bosonizations. 
Thus, it is convenient to use a (simple) continuum de- 
scription of such an integrable model. Interestingly, the 
specific heat behaves as: C^/L = 2ST/{1 + S) + 0{T^), 
where S=l [g[. In the sense of critical theories, such a 
model could be parametrized by a conformal anomaly 
(which is the central charge in the Virasoro algebra): 
C = -j^ = 3/2 p. This fact, together with approximate 



mappings, leads to the conjecture that the criticality of 
this model is governed by operators satisfying a criti- 
cal Wess-Zumino-Witten (WZW) model on the SU(2) 
group at the level k = 2S — 2 |jlffl. Solving numeri- 
cally the Bethe- Ansatz equations the authors of ref. pd| ] 
checked the main conjectures made in refs. [p|,fo|. Thus, 
for (3=1, the effective action for the low-energy excita- 
tions is given by the WZW model on the SU(2) group 
at the level k = 2. The following general statement is 
valid: the Kac-Moody algebra for the group SU(N) at the 
level k=N can be represented by the currents of Majorana 
fermions ^°(a = 1, ..., N'^ — 1) from the adjoint represen- 
tation of this group. Thus, the SU{2)k=2 currents can 
be explicitly rewritten as: I°(I°) = -^'^°''"'^l(r)^l(r)- 
The Hamiltonian is in turn equivalent to the model of 
three massless Majorana fermions ^°(a = 1,2,3) jl^. In 
the continuum limit, the spin-1 chain has only zero sound 
triplet excitations (or magnons, like in conventional mag- 
nets). The staggered magnetization N can be simply ex- 
pressed as N(x) ~ TT\(h + h'^)(T}, where h e SU{2)k=2 
has dimension 3/8 |10| , p^ . 

The spectrum of free conduction electrons yields a sep- 
aration of spin and charge. The spin sector is described 
by an SU(2) WZW model with k=l 0. Here, the two 
fundamental fields g, g^ are spinon pairs which can be 
viewed as free fields except from a purely statistical (in 
that case: semionic) interaction. The electronic spin den- 
sity is usually described as: Sc{x) = Jc(a;)-|-e^*'^^^Nc(a;), 



where J^ = J -f J and N^ ~ Tr{cr(g + g^)} cosv^^c 
are, respectively, the smooth and staggered parts of the 
magnetization. Nc has the dimension 1 and Jc can be ex- 
plicitly written as a function of g and g^ ||7| . The charge 
sector is described in terms of a U(l) scalar field $^ lead- 
ing to a Conformal Field Theory (CFT) with C=l for 
holons as well. Zero sound charge modes propagate with 
the velocity Vp and are ruled by the Luttinger parameter 
Kp = vpjvp. For a ID free electron gas Kp = 1, and 
Vp = 2t sin(fci?a) is the velocity for charge and spin de- 
grees of freedom. The Kondo interaction takes the form: 

H™t=A2(J+J)(I + I) (1) 

-t- aXse'^'' Tr(.gV) cos \/2^$J;N -t- h.c. 
where A2,3 oc Jk, and a is a non- universal constant. 



We have defined i5 = 2kF — tt. Here, we consider sit- 
uations of tlie half-filled case, with (S — > 0. The bare 
conditions are fixed to: Jk <C (J||,t). The coupling 
between the short-range staggering fluctuations becomes 
there prominent. A small term A3 modifies the univer- 
sal dc conductivity of the electron gas CTo = 2e'^KpL/h, 
obtained by applying a static field over a finite part of 
the pure sample. In the presence of magnetic impurities, 
the dc conductivity is rescaled as cr{T) oc {TA3(T)^}~-^ 
[ fi3[ . We assume that the phase coherence is lost after 
each magnetic collision. Then, the diffusion by the array 
of magnetic impurities with S=l provides an anomalous 
power-law, a{T) ex T^/^. Note that the smallest amount 
of Kondo 'disorder' will destroy the so-called Umklapp 
term which gives <7{T) ^ T^^ (when Kp = 1) by using a 
memory- function approximation [|l4| . Clearly, the effects 
of scattering from the array of magnetic impurities be- 
come coherent at lower energy which typically give rise 
to a Kondo localization. An estimate for the localization 
length may be simply obtained from one-loop perturba- 
tion theory; it gives: Cmag = [I/JkT^^. 

The dc conductivity for this state can be easily com- 
puted using the Kubo formula. The resistivity follows 
p{T) ex e'^/'''^^ when ksT < A leading to a perfect 
insulating state. When 5 — 0, the mass gap (which 
can be identified either as the charge gap or the spin 
gap since A3 couples spin and charge sectors) obeys 
A = £~ag = Jk^'^. It is larger than in the usual 
ID KLM (with S=l/2) for the same bare conditions 
p3[. The only possibility to obtain a non-zero conduc- 
tivity is through thermal fluctuations which can acti- 
vate charge excitations. When T ^ 0, correlations of 
cos ^27^$^ show a long-range order. Summing the spin 
part and refermionizing the theory, the diagonalization 
of the charge sector is an easy task even for J 7^ 0. We 
find a gas of spinless fermions ruled by the dispersion 
relation: e(fc) = ±^Jv'^k'^ -I- A^ — vfS/2. It rescales the 
mass gap as A(5) = A — vfS/2 and the insulating state 
can be stabilized as long as A ^ vp6/2. 

Now, let us investigate the spin properties of such an 
insulating Kondo system (for 3 = 0) and explain the 
crucial idea of the paper. The physics behind these cal- 
culations is the formation of bound states between con- 
duction electrons and localized spins, leading to a semi- 
conducting phase. Of course, the occurrence of Kondo 
bound states can be achieved only due to the splitting of 
each S=l localized spin onto two doublet modes. Note 
that although powerful techniques have been used to 
solve the S=l impurity Kondo model, it is difficult to 
have a 'simple' mathematical formulation of the splitting 
phenomenon (see ref. M and references therein). In the 
present lattice model, we express it simply as follows. For 
T ~ A (typically the spin gap) we use the decomposition, 

NAi/« = (Ni + N2) (2) 

with Ni = Tr(/iicr) and hi e SU{2)k=i. Through the fac- 



tor A-'^/^, we consider that the operators hi (with k=l) 
have the scaling dimension 1/2 and not 3/8. Only 'half 
of each S=l spin can form a bound state with a con- 
duction electron. Then, one doublet is assumed to par- 
ticipate in the Kondo process whereas the other would 
remain massless. This phenomenon can be simply de- 
scribed through the important formula: 



^Triga) = A'^/^ Tr(/iIa)Tr(ga) 



(3) 



First, the SU(2) field /i2 remains free (except statisti- 
cal interactions), leading to a critical C=l CFT (i.e. a 
Heisenberg model) with deconfined spinous in the ground 
state. Then, despite of the insulating state the optical 
conductivity would reveal only a pseudogap at A: = tt 
through spin instantons. Following the arguments of 
Tsvelik ||l|], one can define an anomalous fermionic cur- 
rent 3a which is proportional to the topological charge 
density of the staggered magnetization N2. Thus, 3a 
varies linearly with ^e'^'^N2(3^N2 x 3,yN2) which has 
the scaling dimension 5/2. In the (g, w) space, it gives 
a purely dynamical optical conductivity aopt(w, k = tt — 
q) ex — (ui^ — vj^q^)^/^ ~ w^. Second, deconfinement of 
spinons does not imply that conventional magnons fail to 
be quasiparticles in such a spin liquid. The low-energy 
physics also involves optical magnons. It is convenient to 
define: rj = A3(cos V27r$J;(x)) ex A^/^. During the spin-1 
splitting mechanism, the Kondo exchange is rescaled; we 
change, 2ri — + rj. Then, by using Eq.(0), the dimensional 
equality 770 — rjA^^'^ ~ A and omitting high-energy sin- 
glet magnon modes (see later), we find a massive part 
described by: 



Su = 



1 

2^- 



(fx ^TT{VUVU'')-T]o:TiiU + U''): (4) 



where, 7^-^ ex ^vp/ J\\a + y^J\\a/vF, xo = Ty/vFJ\\a, 
xi = X and U{x) = gh\{x) has SU(2) symmetry. 

Excitations of the 0(3) NLcrM are S=l triplets 
(magnons). At low temperatures, Trf7(a;) = constant 
and then magnons become purely massive. Since by con- 
struction TrU has the scaling dimension 1, the spin gap 
(of course) is of order of rja ~ A. The action S(j can be 
rewritten as a triplet of Majorana fields x'^iO' = li2,3) 
with the same mass rrit = A. They mainly contribute 
to a coherent 5-peak to the dynamical susceptibility 
X"(g = 7r-fc >0,uj) ^ mt{q^ + mfy^^"^ 5{uj- y^f~fm^) 
near fc = tt and uj = rrit. On the other hand, the spectral 
density of the staggered magnetization (— l)^N2(x) with 
N2 ^ Tr{h2<j) shows a purely incoherent background. 
Schulz [|l^ explained that it leads to: x"{q > 0;^) — 
sin(2^c)r4-23m[p(f^)p(f^)] with p(x) = ^^^^ 
and c=l/4 for the isotropic Heisenberg chain. Such a 
contribution diverges when uj -^ q a.s (w — q)^^ so that 
the 5-peak due to optical magnons remains distinguish- 
able at u; = i/g^ + mf^. 



During the Kondo localization, the current I is also 
decomposed into two independent parts (Ji,J2) where 
3i are SU(2) currents at the level k=l (with C=l each). 
Some degrees of freedom should account for the miss- 
ing C=l/2. Those are associated with the discrete 
{Z2) nonmagnetic interchange symmetry 1 <-> 2 in the 
definition of I. The marginal term can be rewritten: 
A2(Jc + Jc)(Ji + Ji)- We have the correspondence 
(Ji ® 3c) = ~^<^°''"'XlXl- Then, by using the Dyson 
equation for the fields Xii ^^ fi^^d that the mass rrit is 
sfightly renormalized: A(l-|- ^^^ In ^) and A ~ t. Next, 
we neglect these logarithmic corrections. For J^^a ^ vp, 
free magnon excitations can survive for T < A. We can 
integrate out excitations which correspond to configura- 
tions where TtU 7^ acquires a gap. It becomes an 
0(3) NLcrM where the topological term has no contri- 
bution. Finally, it provides a rescaled mass gap mt ^ 



Aexp(— 27r7 ), where 7 



Jii — t\. Here, the spin 



correlation functions of the massive part are proportional 
to the Mac-Donald function Ko{mtr) with r = (a;o,a;i) 



a.ndKo{x) 



-i/2e-^(l_l/8a;+C'(a;-2)) Q. Mas- 



sive polarons should occur in the interval between nit and 
A due to interactions of electrons and kinks in Ni Q. 

The conclusion that the fixed point is the k=l WZW 
model is consistent with Zamalodchikov's c-theorem [ pT[ . 
This theorem states that if a field theory flows between 
two conformally invariant fixed points, then the value 
of the conformal anomaly parameter C must be smaller 
at the Infra Red stable fixed point. Here, the unstable 
fixed point is ruled by C = 7/2; the stable fixed point 
has C = 1. Deconfined spinous remain as a result of 
the non-Abelian chiral symmetry which is not broken. 
Note that the spinon velocity v '^ J\\a is equal to that in 
the single S=l/2 spin chain problem |l3|. Through the 
transform ^* — > x% magnons become optical. For T <ti 
A, only the gapless modes (spinous) contribute to the 
specific heat Cv = TrT/3v. Now, we choose the particular 
point: J||a = vp- Using the Abelian bosonization we 
proceed to calculate quantities of interest like the spin- 
spin correlations, the NMR relaxation rate 7^ and the 
magnetization by varying the external magnetic field. 

First, notice that an alternative picture of the inte- 
grable S=l chain with /? = 1 is realized in the context of 
spin ladders, when adding a strong four-spin interaction 
l^EjSi(j)Si(j + l)S2(j)S2(j -f 1) equal to the direct 
transverse exchange J± < |1^. Second, using the com- 
mon bosonization scheme for the two-leg ladder (see, for 
instance, Parts II and III of ref. [M), we can rewrite the 
magnetization N of the S=l chain in terms of two scalar 
fields $e and $0, and their duals Og and Qo, as follows: 



Nj, = cos y/nQe cos -s/ttOo 
N.y = sin y/nOe cos \/ttQo 
N^ — sin \/7r*f e cos V^^o 



(5) 



whereas the Hamiltonian of <I>o describes a single free 
Majorana field. It leads to a CFT with a total cen- 
tral charge C=3/2. The operators expiy/7T^e{Qc) and 
cos v^<i>o(9o) acquire the scaling dimensions 1/4 and 
1/8, respectively. Now, we use the representation 

Nc = (cosV2^ec, sinV2^ec, sinV2^*c) (6) 

Then, it is not difficult to show that the splitting of each 
S=l local spin can be still understood through the for- 
mula (||) where Ni=i_2 has the same form as Nc, but 
replacing $c -^ $j and 8c -^ Qi- We have the con- 
straints, $1,2 = l/\/2(<I'e ± $o)- By using the Abelian 
bosonization, Eq.(0) can be understood as follows: only 
one bosonic field, namely <I'i, can hybridize with $c and 
similarly for the dual fields. The term A3 can be simply 
written as: 



'Hint = -771 : cos%/47r$+(a;) : +r]2 : cos V47r<I'-(x) : (7) 
-I- 2ry3 : cos VAttQ^ (x) : 

with the notations: rji ex A/tt^a and $+,- — l/\/2(<I'i ± 
^c). The physics depicted by Hmt becomes similar to 
that of the two- leg ladder system, with V = 0. 

Results obtained in ref. [g^ are then applicable to 
that case and we present them briefly for completness. 
The Abelian bosonization method gives an exact mas- 
sive spectrum which is composed of a triplet of Majorana 
fields x°'{0' = Ij 2, 3) with the same mass mt = A, and a 
singlet x° branch shifted to the high energy nis — — 3A. 
The total N'^+) and relative N*^") staggered magnetiza- 
tion of the massive part are given by. 



N(+) = {aifj,2(73<^, fMia2(T3(7, CriCr2M3Cr} 
N'"' = {MiO'2M3Mi Cri/Z2/Lt3At, Hi 1220-312} 



(8) 



where ai , /i^ are order and disorder parameters of three 
Ising models corresponding to the massive triplet of the 
Majorana fields (x°) while a, jU refer to the singlet branch 
(x°). Masses determine deviations from criticality such 
as mt(^s) = [T — Tc)/Tc >(<) 0. It corresponds to three 
disordered phases with (/ij) 7^ and Uj = 0, and one 
ordered phase with {fi) ^ and (7 = 0. We used the so- 
called Kramers- Wannier duality t —f — t, a ^^ ^ for {(J,^) 
corresponding to a negative mass. Using the asymptotic 
behaviors of the different 2-point Ising correlation func- 
tions pot] and the expression of N2, one obtains: 



(N(-)(r)N(-)(0)) oc ^e-™"' + 0{e 



—3mtr\ 



(9) 



For J± = V, note that $e and Qe are completely free, 



5/2 
(N(+)(r)NW(0)) ~ _Le-2('"*+™=)'- 

(N2(r)N2(0)) ~ - 
r 

For Jk > 0, the S=l magnetization is essentially deter- 
mined by the relative staggered part N^"' = Ni — Nc. 



The singlet branch gives only high-energy corrections so 
that the non-Abelian action Sjj becomes exact at low 
temperatures. The fields $2 and 02 remain free. Thus, 
N2 yields the typical correlation functions of the S=l/2 
isotropic Heisenberg chain. A related quantity of di- 
rect experimental relevance is the NMR relaxation rate 
7i. It is easy to see that the staggered susceptibili- 
ties dominate in 7i~ . Precisely, we obtain: x(2/ci?) = 
1/T + T^^/^ exp — A/T. For a system of non-interacting 
spins xi'^kp) = {TiT)^^ is a constant. Here, the two- 
component underlying magnetic order changes the effec- 
tive field seen by the nuclei and, for instance, T^^ goes 
to a nonzero value as T ^ 0. 

Pleasantly, as in three-coupled spin chain systems [£lj , 
there is a phenomenon which is strikingly analogous to 
the quantum Hall effect. The infinite degeneracy of a 
Heisenberg chain is lifted for a strong magnetic field. 
The system shows a plateau in the magnetization curve 
at magnetization per site m=l/2 and remains insulat- 
ing even when a finite magnetization is induced by the 
external field. Thus, we obtain the same physics as for 
a partially magnetized Stot — 3/2 chain |2^. Regard- 
ing each 5*40* = 3/2 operator as a symmetrized product 
of three S=l/2's, one S=l/2 is polarized by the the ap- 
plied field at each site while the other two form a valence- 
bond-state having a nonlocal topological order associated 
with the breakdown of a hidden discrete Z2 x Z2 sym- 
metry. Here, we have a two-chain realization of such a 
Haldane state and the 'string' order parameter is defined 



as lim|2._j^|_ 



,{Ostr{x,y)) ~ (^1)^(^2)^ 7^ 0. But, we 



expect the vanishing of the Haldane gap for some value 
H = Hci of the magnetic field. We closely follow the 
notations of ref. [p3|. We can see easily that H affects 
only the symmetric field $+ as, Ti+ -^ Ti+ + \/2Hdx^+- 
As long as H < Hci the result is a Sine-Gordon model 
with a gap which is rescaled by the magnetic field. It can 
be viewed as a completely filled lower band for spinless 
fermions with a gap ~ A separating the lower and the up- 
per bands. The Fermi level, lying in the middle, is shifted 
when increasing H. Hd refers to the situation where the 
Fermi energy crosses the gap and just lies at the bottom 
of the upper band. We have the relation A — giisHci- 
Above Hci , the upper band is partially filled leading to an 
amount in the magnetization Sm = m — 1/2. The prob- 
lem is very similar to the commensurate- incommensurate 
transition |24| . Near the plateau 5m oc \J\H — Hci\ and 
1/Ti ~ T^^. The system remains insulating since 0_ 
has a nonzero expectation value and $_ has a gap p3[ . 
In conclusion, we have presented a new Kondo ladder 
system by coupling the so-called Takhtajan-Babujan S=l 
chain to a half-filled ID electron gas and obtained low- 
energy properties using quantum field theory techniques. 
Through the underscreened Kondo effect, this gives novel 
insights into ID spin liquids: the optical conductivity re- 
veals only a pseudogap at fc = tt despite of the insulating 
state, a coherent magnon peak in the dynamical spin 



susceptibility coexists with an incoherent spinon back- 
ground, the magnetization yields a fractional-like quan- 
tization. Such a Kondo state still occurs in the vicinity of 
the /3 = 1 point, but when J\\ (1 — /3) = rut the S=l chain 
is in the Haldane phase |^ and the electron gas is de- 
coupled. These results extend the semi-classical picture 
of ref. ||l|| which is applicable only for Jua/vp —>■ and 
should explain the coexistence between a Kondo singlet 
state and strong antiferromagnetic fluctuations in some 
heavy fermion compounds such as UPt^. More generally, 
this problem belongs to the class of Luttinger liquids in 
active environments, a topic of importance not only to 
the Kondo chain physics, but also to e.g. striped phases 
and carbon nanotubes. K.L.H. thanks B. Coqblin, A. 
Honecker, T.M. Rice and A. Tsvelik for conversations. 
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